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1. Introduction 


Much of economic theory is eurrently diseussed in terms of mathematieal 
eeonomie models. Mathematical economics aims at representation and analy¬ 
sis of problems in eeonomies in order to form meaningful and testable proposi¬ 
tions about eomplex issues often deseribed in a less formal way in everyday life. 
Eeonophysies, on the other hand, originates from attempts at solving problems in 
eeonomies with tools developed by physieists, and is evolving into an interdisci¬ 
plinary research field. Several applieations of the approaeh to stylized models of 
eeonomies have been reeently put forward E,i . The aim of this paper is to show 
how the extremal physieal information (EPI) method of Erieden and Soffer [|^ 0] 
ean be used to devel^ a generalization of the Aoki-Yoshikawa seetoral produe- 
tivity model (AYM) Jlilll]. Below, its modified [|^0] version, which abandons the 


previous, arbitrary metrieal form is presented and the solution to entailed equa¬ 
tions of the fully analytieal formulation of the information prineiples problem [@] 
is given. The approach is based on the maximum-likelihood estimation (MLE) 
and the Eisher information, both the observed and expeeted ones, defined as the 
expeetation value of the observed information widely used in information geome¬ 
try [[sl] and statisties [0,0,1^Si- A similar approaeh was previously used to analyse 
the problem of subjectivity in supply-demand related issues Mm. This paper is 
organized as follows. In Seetion [2] the original formulation of the AYM seetoral 
produetivity model is presented. In Seetion |4] generalization of the method based 
on the EPI method is introdueed. In Seotion[^both approaehes are compared and, 
finally, in Seetion |7]eonelusions are drawn. 


2. The Aoki-Yoshikawa Sectoral Productivity Model 

Sectoral productivity models form key issues in the analysis of productivity 
growth at an intermediate level of aggregation JH]. Such analyses aim to deseribe 
the patterns of produetivity growth aeross and within seetors (e.g. agrieulture, 
manufaeturing, and serviees) and to identify main poliey faetors driving these 
patterns. The natural way of presenting these models is in terms of transition 
probabilities over oeeupation states. In sueh an approaeh, oeeupation veetors and 
partition veetors ean be given an interpretation in terms of eeonomie variables. 
In a mathematical approach, irreducible and aperiodic Markov chains that have a 
unique invariant distribution are used Such an approach makes it possi¬ 

ble to cope with a large number of interaeting heterogeneous agents and to some 
extent ignores the issue of rationality of agents’ behaviour as it is impossible to 
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follow the “motion” of an individual agent in a system eomposed of about 10® 
individuals. Therefore, the assumption that preeise behaviour of eaeh agent is ir¬ 
relevant, is the erueial point. This enables one to adopt some teehniques used in 
statistieal physies and it follows that some models of maeroeoonomies ean be built 
on analogous premises. In their book Masanao Aoki and Hiroshi Yoshikawa 
presented, among others, an interesting model for the eeonomy of a eountry with 
g eeonomie seetors. The ith seetor is eharaeterized by the amount of produetion 
faetor rii, that is, the number of workers in the seetor i, and by its level of produe- 
tivity (effeetiveness) a*. 

Aoki and Yoshikawa were interested in finding the probability distribution of 
the produetivity among seetors. In the statistieal physies language this means that 
the probability distribution of the oeeupation veetor 


n = (ni,n 2 , ...,ng) 


( 1 ) 


of the system is searehed. This eoineides with the standard statistieal physies 
problem of finding alloeations of n partieles to g energy levels. Aeeordirm to 
Boltzmann, the probability distribution of the oeeupation veetor is equal to oD 


n\ 


TT n = 


m=i 


riil 


11 ?"' = 


n\ 


2=1 


Ui=i nil 




( 2 ) 


where p is the probability of the oeeupation of the partieular sth seetor (s = 
1, 2,..., 5 f) by the i\h worker taken to be the same for all partieular eonfigurations 
of these oeeupations. 

The total produetion faetor is exogenously given and fixed so that 

9 

= n , (?7, fixed). (3) 

2=1 


The output of the i-th seetor is given by 


Zi = aini. (4) 

The Gross Domestie Produet (GPD), that is, the total output of the eountry Z is 
equal to 

9 9 

Z = '^Zi = '^aini . (5) 

2=1 2=1 
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In the model it is equal to exogenous aggregated demand D, i.e.: 


Z = D , {D fixed). 


( 6 ) 


Note. It is possible to consider versions of AYM where the demand bound ® 
or the condition of constancy of number of workers ([3]) in the system are relaxed 

iQ. 

The standard Lagrange multipliers method can be used to find the occupa¬ 
tion vector n which maximizes the probability 7r(n) with conserved both the total 
production factor n and the total GPD equal to D. With the help of the Stirling for¬ 
mula ln(f)[^^^ nj) = niilnui — 1) {n » 1 and n* >> 1) the problem 
is reduced to finding the solution of the system of g equations: 


In 7r(f?) -I- v iT'i ~ ~ 

The solution has the following form: 

Hi = n* = , i = 1,2,g . 


A 

drii 


0 . 


(7) 

( 8 ) 


The constants z/ and (3 are determined by inserting ([8]) in ([3]) and ([5])-®. This 
is the Boltzmann distribution for the system which is in the state of the statisti¬ 
cal equilibrium. Scalas and Garibaldi showed iQ that there is a more general 
solution of the form 


rii = n, = 


^—u ppai 


i = 1 , 2 ,...,^; c e 


(9) 


where c is a parameter. Eq.® arises when the appropriate Markovian dynamics 
is taken into account, with the transition probabilities which are tuning, via their 
dependance on the parameter c, the choice of a new productivity sector for workers 
leaving their sector Q. Only in case of c = 0 the Aoki and Yoshikawa solution 
® is recovered and the interpretation of the cases when c 7 ^ 0 can be found in 
0]. 

With the further assumption til]: 


ai = iao i = l,2,...,g, (10) 

where oq denotes the minimal productivity, one gets the most probable vector 
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where r = is the aggregated demand D per agent divided by the minimal 
produetivity. In the limit r 1 one gets 


n,- 




i = 1,2,... ; r > 1 . 


( 12 ) 


The assumption a* = iao also allows for a simplified worker dynamics via cre¬ 
ation and annihilation of components of the occupation vector but due care must 
be taken to make sure that n is conserved 111 ill . Alternatively, in the EPI method 
approach to the AYM model, the probability distribution p (a) of the level a of 
productivity will be found (see Eg. (17^1. If the probability distribution p (f, a), 
which is normalized over space and time fy dt dap {t, a) = 1, then p {t) = 
fy da p (f, a) represents the probability that the worker is found at the time (f, t + 
dt) somewhere within measurement productivity space 3^^, which is a set of pos¬ 
sible values of productivity. Eor example, a high p{t) dt means that there is a high 
chance that a worker is found anywhere in the space of productivities aUhis time. 
In particle physics this property could be called probabilistic creation 


3. Basic information on the information channel capacity 

Suppose that the original random variable Y takes vector values y G 3^ and 
let the k - dimensional vector parameter 9 of the distribution p(y) be the expected 
parameter, i.e., the expectation value of Y: 


9 = E{Y) = [ dyp{y) 
Jy 


(13) 


Eet us now consider the A^-dimensional sample Y = (Yi, F 2 , Yn) = {Yn)n=i, 
where every Yn is the variable Y in the nth population, n = 1,2,N, which 
is characterized by the value of the vector parameter 0^. The specific realization 
of Y takes the form y = (yi,y 2 , ...,y 7 v) = (yn)n=i’ where every datum y„ is 
generated from the distribution Pn(yn|0) of the random variable where the 
d = k X N - dimensional vector parameter 0 [Il2n is given by: 




\N 


(14) 


The set of all possible realizations y of the sample Y forms the sample space B 
of the system. When the variables Yn of the sample Y are independent, then the 
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expected parameter 6*„/ = j^dy P{y\Q) y„/ does not influence the point probabil¬ 
ity distribution Pn(ynl^n) for the sample index n' ^ n. The data are generated in 
agreement with the point probability distributions, which fulfill the condition: 

Pn{yn\Q) = Pn{yn\0n) , wliere n = (15) 

and the likelihood function P{y |0) of the sample y = {yn)n=i is the product: 

N 

P{Q) = P{y\e) = l[pn{yn\0n) . (16) 

n=l 

The Fisher information matrix: Assume that on B the d - dimensional sta¬ 
tistical model: 


S = {Pe^ P{y\e), 0 = e 14 c , 


(17) 


is given, i.e. the family of the probability distributions parameterized by d non- 
random variables which are real-valued and belong to the parametric space 

Ve of the parameter 0, i.e. 0 G 14 C Thus, the logarithm of the likelihood 
function In P ^ 14 —)■ 3^ is defined on the space 14 - 

Let 0 = (0i)f=i ^ be another value of the parameter or a value of the esti¬ 
mator 0 of the parameter 0 = At every point, Pq, the dxd- dimensional 

observed Fisher information (FI) matrix can be defined [ 13. 


iy(0) = -a*4MnP(0) = (-a*4MnP(0) 


le=e 


(18) 


and 0* = djdQi, &■ = djdOi, i, i! = 1,2,.., d. It characterizes the local properties 
of P(j/|0). It is symmetric and in field theory and statistical physics models with 
continuous, regular and normalized distributions, it is positively definite lll3n . We 
restrict the considerations to this case only. T^ expected dxd- dimensional FI 
matrix on S at point Pe is defined as follows 


/^(0)^Pe(lF(0))= / dyP{y\Q) iF{Q) 


(19) 


where the differential element dy = d^y = dyidy 2 ---dyN- The subscript 0 in the 
expected value signifies the true value of the parameter under which the data y ar e 
generated. The FI matrix defines on S the Riemannian Rao-Fisher metric [fs. 121. 
Sometimes, due to the probability distribution normalization and the regularity 
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condition, the dx d - dimensional observed Fisher information (FI) matrix ean be 
recorded in the "quadratie" form [[sl]: 

IF = (^a*'lnP(0) aMnP(0)) . (20) 

The eentral quantity of EPI analysis is the information ehannel eapaeity / whieh 
is the traee of the (expeeted) Fisher information matrix. Sinee under above eondi- 
tions, the observed Fisher information matrix is diagonal iF(0) = diag(iF„„(0)), 
henee the information ehannel eapaeity /(0) is equal to: 


N 

/(e) = V / <;j//>(!,|e)iF„„(e) 



( 21 ) 


where i 


P(0) J2n=i ^nn(0) Is the information channel density HH. 


4. Generalization of the Aoki-Yoshikawa Model 


The generalization of the AYM presented in this paper eonsists in eonsidering 
productivity as a continuous random variable A. The transition from the discrete 
variable to the continuous one is performed via 

72 ■ 

A = tti ^ a and Pi = — ^ p(a) . (22) 

n 

As a consequenee, the probability distribution function p has to be normalized: 


9 . 

'y^^Pi — dap{a) = 1 , 


(23) 


where fa denotes the set of possible values of the produetivity. Analogously, the 
expeetation value of the produetivity is replaeed in the following way 

9 . 

Oa = (A) = yy^pi Oj ^ / da p{a) a , (24) 


with the eonstraint 


(A) = D/n . 


(25) 
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In order to find the probability distribution of the level of produetivity A the 
EPI method is used. Below, the forms of the Fisher information will be adopted 
to estimate the sealar parameter 0a, (l24l) . lll4l] : 


0A = (A) = da p{a) a , 


(26) 


’ya 


whieh in this paper is the expeetation value of the random variable of the pro¬ 
duetivity level A. The basie information on the Fisher information is given in 
Seetion_p and the eonstruetion of the information ehannel eapaeity ean be found 
in&M. 

Previous attempt to solve the Aoki-Yoshikawa produetivity problem by the 
EPI method was based on the Frieden-Soffer approaeh J^]. However, this paper 
implements only the analytieal form of the (modified) observed struetural infor¬ 
mation prineiple [Il4l] . 


4.1. The expected Fisher information and the information capacity of the chan¬ 
nel {6 a) 

Aeeording to the the main assumption of the EPI method proposed by Frieden 
and Soffer, the system itself samples the spaee of positions (that is, the spaee 
of values of produetivity A levels). This spaee is aeeessible using its Fisherian 
degrees of freedom [H, 12]. The sample of the EPI method (so-ealled inner sample 
Hi , 14]) for the AY model is = 1 -dimensional. Thus, both the sample spaee 
B and the base spaee of events are in AY ease equivalent to ya H. 


11,111] the condition of the minimal value of the information (kinematical) channel capacity 
/ —min is postulated as the one that hxes the value of A in a unique way. However, sometimes 
the non-minimal values of I are also discussed as they lead to the EPI method’s models, which 
are of a physical significance Si. Some discussion on this topic can be also found in iflj . For 
example, from the analogy with the EPI model of momentum distribution presented in i, for 
A > 1 the nonequilibrium, stationary solutions for the square flow (where the flow is proportional 
to the root productivity) might be obtained instead. The point is that when choosing the size A 
of a sample, different classes of the EPI models are obtained. Yet, with the value of A fixed, 
the particular EPI model is provided by the solution of one or two information principles. In the 
case of the discussed AY model, the observed structural principle, which is a differential one, 
and the variational principle are used. The variational information principle is connected with 
the extremization of the physical information (see Section IT2l l. Thus, the condition I —?> min 
is choosing the type of the model only IJ, [tI [iS]] whereas with the information principles, the 
particular solution inside this type of the model is found. 
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As = 1 thus p{a\6A) = p{o) is the likelihood function for the AY model and 
the d = 1-dimensional statistieal spaee of for the AY model is as follows: 


= {P0A = I^a), Oa e Ve^ C 


(27) 


where is the parameter spaee of 6a- For N = 1 the information ehannel 
eapaeity I reduees to the Fisher information IpiOi) = If{0a) for 9i = 9 a, the 
only parameter. The neeessary steps leading from the general form of the Fisher 
information to the one used in this paper are similar as in H. 

The probability amplitude Qa = q (a|6*A) is defined in the following way il6 . 

iH: 


p{a\^A) = 4^^ («I^a) 


(28) 


The N = 1-dimensionality of the sample means that the rank of the amplitude 
q {a\9A) of the (produetivity) field is also equal to 1 lldl]. 

As 9a is the sealar parameter and the dimension of the sample is equal to 
= 1 the information ehannel eapaeity I{9a) of the measurement channel (9a) 
112] is equal to the Fisher information Ip(9a) of the parameter ^aIQ: 


1(9a) = Ip(9a) . 


(29) 


If((^a) is the information about the unknown parameter 9a eonfined in the A^ = 1- 
dimensional sample for the random variable A. 

The analytical form of the (expeeted) Fisher information on parameter 9a is 
equal to [14]: 


If(0a) = / dap(a\9A) i\Fa(9A) 


'ya 


da p (a \9a) 


'ya 


hip (a \9a) 


d9\ 




'ya 


d9\ 


p(a\9A) V 99a 


'ya 


'ya 


d^p(al9A) , k 

„ d^p(al9A) 

da ( - qaqa + -- 


dp (a\9 a) 


(30) 
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where Eq.(l28l) and the denotations q'^ = and g” = have been 

used and the index in iFa(6*^) signifies the dependence of the observed Fisher 
information on a. In the last equality the relation: 


d'^p{a IOa) 

del 


1 

2 


{Qa? + ^(laqa 


(31) 


was also applied. 

Due to the normalization, (l23l) : 


dap{a\9A) = 1 , 


(32) 


'Va 


and the regularity condition OQ: 

^ d'^p{a\9A) 0^ / 


da 


'ya 


de\ 


Jy. 


dap{a\9A) = -^1 = 0 , 
o9a 


(33) 


the analytical form (l30l) of the Fisher information transforms into the following 
metric form 11411 : 


If{9a) = / dap{a\9A)iiFa{9A) 

Jya 


da 


ya P{a\0A) 


dp {a \9a ) 

d9. 




(34) 


Due to Fq. (l3^ and in accordance with Fq. (l29l) . both the EPl method form of 
the (expected) Fisher information [Il4l] for the AYM and its information channel 
capacity for the measurement channel (9a) are equal to: 


I(9a) = If(9a) = -[ da qa{9A)ql{9A) ■ (35) 

Jya 

The information channel capacity I is the one which enters into the estima¬ 
tion procedure of the FPI method. The presented below derivation of the form 
of amplitude qa as the self-consistent solution of the information principles con¬ 
sistently uses the analytical form iHH of the structural information principle, 
and in this respect, it is different from the derivation of the Boltzmann distribution 
given in |3]. 


10 











4.2. The information principles and generating equation 
In [13] the existence of the (total) physical informatio^ K 

K = I + Q>{) (36) 

was postulatec0. The choice of the intuitive condition ii' > 0 is connected with 
the expected structural information principle of the EPI method: 


/ + K Q = 0 


(37) 


derived for k = 1 in [[60, where k is the so-called efficiency coefficient introduced 
in [3]. The general forms of / and the structural information Q are given in [3,0] • 
The form of the information principle, which is more fundamental than ([3^. is the 
observed structural information principle that has the form cf + iF = 0 Jg] . The 
derivation of the particular form of the observed structural information principle 
for the AY model is similar to the one for the EPR-Bohm problem ItldH so only 
necessary steps will be presented. 

The other information principle of the EPI method is the variational informa¬ 
tion principle [3]: 


5{i + g) = 0 


(38) 


It has to be stressed that it is the (modified) observed structural information prin¬ 
ciple fli, HQ, and not the expected one, which is solved self-consistently 
together with the variational information principle. Below, both information prin¬ 
ciples will be constructed and solved in case of the AYM. 

Note. To obtain the value of the efficiency coefficient n, the information 
principles, i.e., the (modified) observed structural information principle HQ 
and the variational information principle, together with the physical preconditions 
which are specific for the model, e.g., some symmetry conditions, have to be 
solved simultaneously [14 h As a result both, the specific form of Q and the value 
of K are obtained il. Inli it is suggested that for the EPI models which have the 
quantum counterparts, n = 1 IQ . whereas for the classical models 0 < k < 1 
[ 3 ]. Prom the below analysis it follows that in the Aoki-Yoshikawa model k = 1, 


^With such general understanding of K, the diversity of the equations of the EPI method is a 
consequence of dive rse p reconditions dictated by the investigated phenomenon ll3ifl[l^l^). 

®See also iSlIll 14]l . where the differences between the Frieden-Soffer original form of the 
physical information and information principles used in this paper are discussed. 
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analogously as in the EPI model of the Boltzmann energy distribution |3] • From 
the analysis presented in [Il4ll it also follows that Frieden’s EPI method of “eov- 
erage of quantum meehanies” ean be eonstrueted by giving the quantum meehan- 
ieal interpretation to the statistieal probability amplitudes. Yet, in the ease of the 
Einstein-Podolsky-Rosen-Bohm problem the quantum eharaeter of the amplitudes 
is provided by the EPI statistieal information theory modelling itself H. Thus, 
the amplitudes QaidA) of the EPI method for the Aoki-Yoshikawa model eould 
gain the quantum meehanieal interpretation if only a reason for the quantization 
(e.g., as tti = iao in Eq.(fT0l)l of the produetivity levels exists, e.^ as in the Ein¬ 
stein model for speeifie heat lilSll . The assumption (fTOl) is used in yj] for the AYM, 
whieh is quoted only in order to eompare this model with the one obtained via the 
EPI method in Seotion[^ 


4.2.1. The information principles for the Aoki-Yoshikawa model 

The struetural information Q in the AYM for the system deseribed by the 
set of amplitudes [Il4ll is as follows: 




da ql{9A) 


(39) 


'Va 


where for the simplieity reason the denotation = qa{dA) = g(a| 6 *A) 4 S^sed. 
Now, the physieal information K, (l3^ . in the AYM is as follows 


K = I + Q= / da kaiOA) , 
Jya 


(40) 


where I is given by Eq.(l30l). In Eq.d^. ka{9A) is the density of the physical 
information, whieh aeeording to Eqs. (l30l) . (l39l) and (ISTl) is equal to 


, //i \ " (^^p(a|6'A) 1 2 / \ 

kaiOA) = - qaQa + - ^^aiQa) 

= -\(laqa + ^{qy + ^ql^a{qa) 

= —qaq'a + ^ql'^aiqa) , 


(41) 


where the modified observed structural information used in the EPI method 
has been introdueed [Il4ll : 


Cf’a(<?a) : = 


qli^A) 


{qaf + ^a{qa) • 


(42) 
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Under the assumption of analytieity of the log-likelihood function hip (a |6 *a ), 

6 a I around the true value of 9 a with the use of the 


the Taylor expansion Inp ya 

denotations and 

following form of the observed structural information lll4l] : 


, leads to the 


'dA=6A 


^aiQa) = 


1 


&a) 


2 qa{9A)q:{eA) - {q:{9A)y 


(43) 


Here the appearance of qa in the argument of means that the probability 
p{a\9A) (and its derivatives) present in in the derivation [|T3l of (1431) . has been 
replaced by the amplitude Qa (and its derivatives). 

In what follows, the forms of the amplitudes Qa that are the solution to the 
AYM will be searched for among combinations of the exponential functions. Ad¬ 
ditional assumption that the term with the first derivative q^^id a) on the RHS of 
the above equation cancels with a term in ^a{qa) has been made. 

Now, after moving the term f in Eq. (l43l) from the Fisher information 
part on its RHS to the structural one on its LHS the modified observed structural 
information principle was obtained [Il4n . (This shift between c^aiqa) and c^aiqa) 
is then used in Eq. (l4T]) .l Thus, the modified observed structural information prin¬ 
ciple for the AYM has the following form lll4l] : 


- 2 qa{6A)qa{6A) + ql{dA) c^aiqa) = 0 , (44) 

where 

^a(ga) = (^C^aiqa) + -^^^^‘2{qa{6A)f'^ 

= ^7^2g„(0A)gl'(0A) • (45) 

qa[6A) 

Equation (l44l) arises purely as a result of analytieity of the log-likelihood function. 

The EHS of Eq. (l44l) is (up to the factor 3) the density of the physical informa¬ 
tion ka{d a) given by Eq. (l4TI) . This one is the function of the observed structural 
information c^aiqa) (which at most can be the function of the amplitudes qa{6A)), 
of the amplitudes themselves qa{6 a) and of their second derivatives. 

Comment: In the AYM, the efficiency factor n is equal to k = 1 |3]- This fol¬ 
lows from the fact that except for the information principles, no additional differ¬ 
ential constraints are put upon the amplitudes Thus, the presented EPI model is 
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a pure analytic one |@], similarly in this respect as in the EPR-Bohm problem [[Till . 
Using Eq. (l4TI) the physical information K, (|4^ . takes the following form: 

K = I + Q 

= da 

JVa 

Erom Eq.dUl) the expected structural information principle (see Eq. (IT7]) l. for k = 
1 , follows: 


(46) 


\ \ 

-^fladi + ^ql^OA) ^aida) 


/ + Q = 0 , (47) 

where / + Q is given by the RHS of Eq.(l46l). 

The differential equation (l44l) is the first one from the information principles 
used in the EPI method. The second one presented below is the variational infor¬ 
mation principle 

In order to obtain the variational information principle, we have to transform 
the physical information K, (l46l) . into the metric form, i.e., the one quadratic in 
q^. Therefore, after integration by parts, K can be rewritten as follows iQ: 

K = I + Q= f da(^kr\0A)-Y) ^ (48) 

where the constant Ca is equal to: 

Ca = (qa (oo) q'^ (CX)) - qa (oo) q'a (ao)) , (49) 

where Oq is the smallest (absolute) level of the productivity and is the 

metric form of density of the physical information: 

C"(»a) -fk + l AVa) ¥,(9a) . (50) 

The variational information principle has the form lilii: 

= 6(,.)(I + Q)= (51) 

= = 0 . 
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The solution of the variational problem (|5T]) with respect to Qa is the Euler- 
Lagrange equation: 


d _ dk^^\eA) 

dOA V dq'ai^A) J dqa 


From this equation and for k'^^^{9A) as in Eq. dSOl) . the following differential equa¬ 
tion is obtained for every amplitude Qa'. 


“ 2 dq^ 


(53) 


As ql{0A)'^a{c[a) ts explicitly the function of qa only, the total derivative has 
replaced the partial derivative over qa present in EqJ52l). The obtained form of 
equation (15^ differs slightly from the Erieden form and is the same as in [Il4ll . 
The origin of this difference is the fully analytical form of density of the physical 
information (|4T]) . 

The modified observed structural information principle (l4^ and the varia¬ 
tional information principle (ISTI) (from which the Euler-Eagrange equation (15^ 
follows) serve for the derivation of the equation which generates the distribution. 


4.2.2. The derivation of the generating equation 
Using the relation (l5^ in Eq.(l4^. one can obtain: 

1 d{q^^a{.qa)) 2 


o 'a 
2 dqa 


= Qa^aiqa 


The above equation can be rewritten in a handier form: 

2dqa _ d (Iql^aida)) 


qa 


ln2 


ql<^a{qa) 


(54) 


(55) 


from which, after integration on both sides, the following results can be obtained: 


^ql{dA)<^a{qa) = C? ql[9A) 

hence okqa) = 2a^ , (56) 

where the constant of integration c? is a complex number in general. By sub¬ 
stituting Eq.(l5^ into Eq.(l5^. we obtain the searched for differential generating 
equation for the amplitudes qa [0]: 

q'aiOA) = qai^A) for Qa^Vq^, (57) 
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which is the consequence of both information princMes - the structural and vari¬ 
ational ones. This result was obtained previously in [[^ for the Boltzmann proba¬ 
bility distribution but the arrival at the structural information principle is here H 
different and the form of both information principles also differs slightly. 

Note. If an explicit dependence of c^aila) on the productivity a is assumed, 
i.e. 0,), then a wider scope of solutions to the problem (l54l) is possible, 

which also includes non-equilibrium solutions 10] ■ These solutions correspond 
to the non-equal probability of the occupation of the particular sth sector by the 
ith worker for all particular configurations of these occupations (contrary to the 
assumption used in (|2l)). 


5. The probability distribution for the Aoki-Yoshikawa model 


5.1. The definition of the variable of the additive fluctuations 

The EPI method analysis for the distribution of the level of productivity is 
in accord with the general approach of Frieden. The displacement Xa, defined 
as Xa = A — (A), is used instead of values of the productivity level A. Thus, 
the additive partition is performed: Va = A = (A) -f Xa (similarly, as for the 
Boltzmann distribution in 11). It can be performed at the level of the information 
channel capacity, as it was originally proposed in H and developed in [12] for 
the general distribution which is free of necessity to set the requirement for the 
shift-invariance, or it can be made at the level of the generating equation. The 
latter possibility has been chosen as in the considered case this is the simple one, 
i.e.: 


ya = a = 0A A Xa , ao <ya <oo , = ao - hA < Xa < oo , (58) 

where Xa = is a particular displacement. The simplifying assumption that the 
fluctuation of productivity is unbounded from above, is used (compare [H for the 
discussion on distribution of the energy fluctuation). Then, the EPI model is built 
over the space Xa of the displacements x^, which in our case is [flill . 

A simplifying notation now will be introduced: 

qe^ixa) = q(xa + 9 a\0a) = q{a\dA) , (59) 

which leaves the whole information on 6 a that characterizes q{xa + 6A\dA) in 
the index of the amplitude qg^{xa) (and similarly for the original distribution 

PbaM = p{xa + 6a\6a) = p{a\9A)). 
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Now, appealing to the “ehain rule” for the derivative: 


d d{a — 9 a) 


d 


d 


dOt 


dOA d{a — 6 a) d{a — 6 a) 


cix„ 


( 60 ) 


a transfer from the statistieal form (157]) of the generating equation to its kinemati- 
calforii^ 


d^qg^{:)^a) 

d^l 


a'^qe^i^a) , 


(63) 


where qg^{^Xa) is the amplitude of the distribution of the produetivity level fluetu- 
ation and it was ehosen a to be a real eonstant (see footnote [8]). 


5.2. The solution of the generating equation 

As the amplitude qg^^ is a real one, thus in Eq.(l63l) has also to be real. When 
the value of the fluetuation of the produetivity is not bounded from above, and 
this eondition is realized by approaehing infinity, then a has to be real. In 
this ease from the normalization of the squared amplitude we get 


1 

4 



cixaPe^(x„) = 1 , 


(64) 


it follows that the solution of Eq.(l63]) is purely of an exponential eharaeteiFI fll : 


qg^ (Xa) = B exp {-a x^) + C exp (a x^) , a G R+ , (65) 


^ Note; Taking into account that dxa = dya, which is connected with the fact that parameter 
9a is a constant, we can transfer from the statistical form of the physical information K = I + Q 
(l40l) - (l4T]) to its kinematical form with the information channel capacity as follows 


I = 



/ dqg^iKg) 

V dxa 


2 


and the structural information in the form: 



(Xa)'^ Cfg{qg^ (X„ 


)) 


( 61 ) 


(62) 


®The other possibilityfor a is that it is a purely imaginary number. Then the solution has the 
trigonometric character 1^. Yet, when oo then due to the normalization condition (l64t . 

the trigonometric solution is not the allowed one. (The case when the parametric space is finite 


can lead to the trigonometric solution, as it is in case of the EPR-Bohm problem Il4l ll4ll l. 
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where B and C are real constants. As the normalization condition (l64l) is defined 
in the interval < cx) thus, the part of the solution with the positive 

exponent has to be rejected due to its divergence to infinity. Therefore, the re¬ 
quirement that a > 0 leads to C = 0. 

In summary, the searched form of the amplitude is as follows: 

qeA^a) = B exp (-ax„) , a e R+ , x™” < x^ < oo . (66) 

From this and from the normalization condition (l64l) . the constant B is obtained: 

B = ±2 \/^ exp (ax™’*) . (67) 


Thus, the final form of the amplitude for a G R+ has the form 

70 a (Xa) = ± 2 exp [a (x™” - x^)] , a G R+ (68) 

and a is given in units [1/productivity]. 

The following probability distribution of the fluctuation of productivity level 
Xa can be obtained from the amplitude (1681) : 


p (x„) = (Xa) = 2q! exp [2q! (x™" - Xa)] , a G R^ 


(69) 


In accordance with Eq. (l5^ the following relation holds a = 9 a + thus, 
da/dy:.a = 1. Therefore, the distribution of the random variable A has the form: 


p{a)=p (x„) 


I da / dx, 

Now, as the expectation value of A is equal tc@: 


= 2a exp [—2a (a — oq)] , qq < a < oo . (70) 


+ 00 


(A) = 9a = dap{a) a , 


ao 


thus, inserting Eq.dTOl) into (ITT]) , we obtain: 


(71) 


2a = {{A) - ao) 


-1 


(72) 


® Let us notice that from Ea. dTTl i it follows that A is the unbiased estimator of the expectation 
value (A) of the level of productivity, i.e., (A) = A. 
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Figure 1: This figure depicts the comparison of the cumulative probability distributions 
= p (a) da calculated for the probability distribution of the level of productivity given by 
Eg.dTTli with the one observed across workers (W) [source ifl^ for the 2005 year’s data] (thick 
dotted line). The productivity cut “a” is given in the unit of 10® yen/person in and on both 
axis the common logarithmic scale is used. The smallest absolute level of the productivity oq of 
the worker is taken to be equal to zero. The ratio D jn (which value is equal to the Boltzmann 
temperature ifT^ ) has its central value equal to 135 (long dashed line). For comparison, two 
additional curves, the left one with D/n = 100 (solid line) and the right one with D/n = 170 
(short dashed line) are also plotted 01911 . 


Taking into account the constraint (^4) = D/n, (l25l) . using Eq. d^ finally the 
amplitude ean be obtained: 


(IdAK^a) 


= ±2 




exp 


n — ao 


Xg + {D/n - gp) 
2{D/n — ao) 


(73) 


where ao — D/n < Xg < oo 

and the seareh for probability distribution of the level of produetivity: 


p (a) = (D/n)-ao 


— exp 

—n.r, ^ 


CL — (iQ 

{D/n)-ao 


for 

for 


a> ao 
a < ao 


(74) 
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The distribution (1741) is the final result of the EPI method for AY model of pro¬ 
ductivity. 

Finally, Figure 1 shows the comparison of the cumulative probability distri¬ 
butions P^\a) = p {a) da considered in [19] for the probability distribution 


of the level of productivity given by Ec 


741) with the observed productivity dis¬ 


tribution across workers (W) (source P19n i. If one considers the whole range of 
the productivity cut “a” then the exponential law (1741) (with the smallest absolute 
level of the productivity oq of the worker equal to zero) fits the data reasonably 
well. 


6. Comparison of two approaches 

To compare the approaches, the discrete variable and assumptions have to be 
reproduced. The assumptions that have been made in the AYM are the following 

dH): 


ai = i ao where i = 1, 2,..., , 


(75) 


the number of production sectors is great, >> 1, and 

D/n 
Oo 


(76) 


This leads to: 


F(l|n*) 



* = 1,2,...; r»l,(77) 


where n*, i = are the coordinates of the most probable occupation 

vector n dU). The above equation gives the probability that a randomly selected 
worker is in the ith sector, provided that the economy is in the state n* = (n*, n^, 



Therefore, to compare both methods the integration of the obtained probability 
distribution (1741) in the segments (oj, Oj+i) is indispensable. Firstly, let us consider 
the case for which the "width" of the sectors is equal to oq = the smallest 
(absolute) level of the productivity. With this assumption a* = * amin, where 
Omm = ao > 0, (IT^ . the result is obtained: 


»(i+l)ao 


Pit) = 


da p (a) = (1 - for i = 1, 2,... , (78) 


' laQ 
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which in the limit r >> 1 gives: 


P [i) PS (1 + ^—) I e ^ ) for i = 1,2,... and Oq > 0 , r » 1 .(79) 

r 2r^ \ r J 

Seeondly, let us eonsider the ease of oq = 0 with the "width" of the seetors equal 
to 5a. This leads to 

pi5a 

P{i)= da p (a) = (—1 + i = 1,2,... for oq = 0 , (80) 

J {i—l)5a 


where 


D/n 

r = 


5a ’ 


(81) 


has been introdueed instead of r (Om . In this ease for r >> 1 one ean finally get: 
P ii) - (^ + f = 1, 2,... , for ao = 0 ; f >> 1 . (82) 


This means that for large f both methods give the same results in the first order 
approximation. Note that the solution (l80l) is exaet and in the AYM approaeh 
some additional assumptions have to be made to obtain the final formula. 


7. Conclusions 

The Aoki-Yoshikawa model, although relatively simple, gives interesting re¬ 
sults and ean be used as a starting point for various analyses. Here we have 
adopted methods used in theoretieal physies to generalize the model and as sueh 
these methods refer to the same phenomenon as the original one JH] (see also 
Figure 1 in Seetion [5^ and the diseussion in [lllll l. Our approaeh allows for ex¬ 
aet solutions. The original AY model and the approaeh presented in this paper 
agree in first order approximation. Models of phenomena eonstrueted within the 
proposed approaeh ean be used as tests of the Extreme Physieal Information Prin- 
eiple and we envisage sueeessful applieation of these methods in other fields of 
researeh [fifllil. 

The original EPI method was invented by Frieden and Soffer [1^3]. They, to¬ 
gether with Plastino and Plastino, put the solution of the (differential) information 
prineiples for various EPI models into praetiee Nevertheless, the derivation of 
the generating equation (1571) for probability distribution of the level of produetivity 
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(r74b (which could have been inferred by the eomparison with the Boltzmann dis¬ 
tribution [0]) differs from the one used in the original Frieden-Soffer approaeh [0]. 
The main differenee in the presented derivation is that in this paper the observed 
physieal information used direetly in the struetural information prineiple was eon- 
sistently obtained from the analytieity eondition of the log-likelihood funetion 
0,0] without any jump from its analytic to its metric form. Only then, the gen¬ 
erating equation (1571) and (1^ was derived. This allows eonseeutively to obtain 
the probability distribution (1741) of the level of produetivity for the statistieal in¬ 
formational generalization of the Aoki-Yoshikawa seetoral productivity model. 
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